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Abstract 

In the preceding paper [1] we dealt with discrete signals of finite duration. Here we 
generalize previous results and demonstrate that the Dirac representation theory can be 
effectively adjusted and applied to continuous or discrete signals of infinite time 
duration. The role of the identity and projection operators is emphasized. The sampling 
theorem is viewed from the point of view of orthogonal physical states. An orthogonal 
basis which spanned the time space, ceases to be orthogonal and becomes overcomplete 
if the domain of frequencies is restricted in a bandwidth. In this case there exists an 
infinite number of sub-bases of discrete times which are orthogonal and complete. The 
relation between the overcomplete bases and a complete one is the essence of the 
sampling theorem. The signal theory is reformulated in the framework of the Dirac 
bra-kets. The case of signals existing for positive time is treated in detail. 
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1. Introduction 

The usual approach to signal theory is to start with the continuous case and treat the 
discrete finite one as an approximation. In the preceding paper [1], which we will 
denote by [I], we have adopted a different approach for the following reasons: 

1. In practice the data are collected in finite discrete sequences [2]. 

2. The discrete finite case has its own peculiarities, different from the 
continuous case. Moreover, there are several, possibly many, discrete cases for 
each continuous case. The continuous case, if it exists, is a common limit of 
different discrete cases [3-5]. 

3. The transition to the continuos case is more complicated, as the limits can be 
distributions (generalized functions) [6-8]. 

Most of the developments in this paper are based on the assumption of Dirac [9], 

made in 1930, that a given self-adjoint operator A can be presented as 

A= j a|a)(a|da, (1.1) 

a(A) 

where A| a) = a| a) , and a(A) is the spectrum of A . Dirac assumed that all properties 

of hermitean operators on finite dimensional vector spaces can be extended and be 
valid in infinite dimensional Hilbert spaces even for unbounded self-adjoint operators. 
A full mathematical justification of the Dirac formalism came only about thirty years 
later. At the beginning the formalism was severely criticized by mathematicians. 

Von Neumann in the preface to his 1932 fundamental book [3]: "Mathematical 
Foundations of Quantum Mechanics" wrote: "Dirac, in several papers, as well as in his 
recently published book, has given a representation of quantum mechanics which is 
scarcely to be surpassed in brevity and elegance, and which is at the same time of 
invariant character. It is therefore perhaps fitting to advance a few arguments on behalf 
of our method, which deviates considerably from that of Dirac. 

The method of Dirac, mentioned above, (and this is overlooked today in a great 
part of quantum mechanical literature, because of the clarity and elegance of the 
theory) in no way satisfies the requirements of mathematical rigor ~ not even if 
these are reduced in a natural and proper fashion to the extent common elsewhere in 
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theoretical physics. For example, the method adheres to the fiction that each 
self-adjoint operator can be put in diagonal form. In the case of those operators for 
which this is not actually the case, this requires the introduction of "improper" 
functions with self-contradictory properties. The insertion of such a mathematical 
"fiction" is frequently necessary in Dirac's approach, even though the problem at 
hand is merely one of calculating numerically the result of a clearly defined 
experiment. There would be no objection here if these concepts, which cannot be 
incorporated into the present day framework of analysis, were intrinsically necessary 
for the physical theory. Thus, as Newtonian mechanics first brought about the 
development of the infinitesimal calculus, which, in its original form, was 
undoubtedly not self consistent, so quantum mechanics might suggest a new 
structure for our "analysis of infinitely many variables" ~ i.e., the mathematical 
technique would have to be changed, and not the physical theory. But this is by no 
means the case. It should rather be pointed out that the quantum mechanical 
'Transformation theory" can be established in a manner which is just as clear and 
unified, but which is also without mathematical objections. It should be emphasized 
that the correct structure need not consist in a mathematical refinement and 
explanation of the Dirac method, but rather that it requires a procedure differing 
from the very beginning, namely, the reliance on the Hilbert theory of operators." 

The above citation was written before the invention of distributions [lc],[ld] in 
the 1950's. The theory of distributions gave new insight into the Dirac formalism 
and relieved it from the accusation of using improper functions. A complete 
justification of the Dirac formalism was given by Gel'fand [10-12], who introduced 
the rigged Hilbert space and the generalized eigenvectors [10-13]. Within the 
framework of this theory, eq. (1.1) is fully justified and can be applied in a simple 
manner, which generalises the finite-dimensional results of paper [I]. The theory of 
Gel'fand is also the basis of axiomatic quantum field theory [14] and was further 
adapted and elaborated in order to describe more plainly quantum mechanics in the 
Dirac formalism [15-20]. 

In the treatment of the finite discrete signals we used the bra-ket formalism [9] and 
ket bases to describe the transition from signals represented by time points to their 
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frequency representation. As a guideline we used the orthonormaHty condition, which is 
the condition that physical states (eigenstates of self-adjoint operators) satisfy in 
quantum mechanics (in one dimension, and for different eigenvalues). The peculiarity 
of using ket-vectors and not wavefunctions is that orthogonality can be employed (at 
least for the problems of this paper), without the need to explicitly refer to boundary 
conditions. One of the reasons is that wavefunctions in Hilbert space are required to 
vanish at space infinity, while for ket-vectors (in the rigged Hilbert space) space infinity 
is a regular point. We continue here the same approach and consider, for example, the 
sampling theorem from this point of view, and we find that the use of orthogonality and 
projection operators lead us to understand the sampling theorem as coming from a 
relation between an overcomplete non-orthogonal basis and an orthonormal complete 
one. Other aspects of signal theory are treated as well in the bra-ket formalism. 

In this paper we introduce the Dirac representation as a framework or a basis for 
signal theory. As the subject is quite broad we will concentrate in the paper only on 
theoretical considerations barely touching applications, leaving them for next 
publications. 

2 The infinite dimensional Dirac space 

In paper [I] we considered a basis with a discrete and finite index. The treatment of the 
basis with a continuous index is more involved. For the continuous case the ket basis 
(13.4) can be generalized to a ket basis labeled by the time t : 

|t), where: -oo<t<oo, (2.1a) 

with the dual bra space: 

(t|, where: -oo<t<oo, (2.1b) 
for which there exists the scalar product (in terms of distributions [6-8]): 

(t|t') = 5(t-t'), (2.2) 
where 5 is the Dirac delta function (see Appendix A). The identity operator now takes 
the following form: 



i= j |t)(t|dt. (2.3) 

— oo 

Equations (2.1-3) define the largest Dirac space for continuous signals. In a way similar 
to eq. (1-2. 1), a ket- vector and a bra- vector in this space (a rigged Hilbert space ) will 
have the most general expansion with the bases (2. la) and (2. lb) in the form: 

oo oo 

|u)= j u(t)|t)dt, (u|= j (u(t))*(t|dt, (2.4) 

— oo — oo 

respectively, where u(t) are the expansion coefficients, which can be distributions. From 
the orthonormality condition (2.3) one gets: 

(t|u) = u(t), (u|t) = (u(t))*. (2.5) 

Let us take as an example the signal f(t) , which can be considered as a ket | f ) with 

projections into the orthonormal basis |t) (i.e. the components) equal to the values of 

the function, i.e.: 

(t|f) = f(t). (2.6) 
In the Dirac space one can define other orthonormal ket-bases, among them the one 
labeled with the angular frequency co: | co) . Its transition to the t-basis is obtained via the 
matrix elements similar to eq. (1-3.8): 

{t|co) = ^e'", (2.7) 

from which we can get that 

I«> = 7=1 e'"|t>dt. (2.8) 
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The normalization in eq. (2.7) was so chosen to ensure the orthonormality condition 

(co|co') = — f e'("-'"')'dt = 5(co-co') . (2.9) 

— oo 

As the co-basis spans the same Dirac space as the t-basis, the identity operator satisfies: 

oo oo 

i= j |t)(t|dt= j |co)(co|d(0. (2.10) 

— oo — oo 

Similar to eq. (1-4.6) we can introduce the angular frequency operator: 



00= j co|co)(co|dco, (2.11) 

— oo 

having the following properties 

cb|co) = co|co), (co|cb = (co|co, cb"|co) = co"|co). (2.12) 
Using eqs. (2.10) and (2.7-9) one can prove that 

oo oo 

(t|(ia))|u)= j (t|(ia))|co)(co|u)dco= j ico(t|co)(co|u)dco 



If d 7 d 

= ^= icoe'™(co|u)dco = — (t|co)(co|u)dco = — (tlu 



(2.13) 



i.e. in that sense the ico operator can be considered as the time derivative operator. 
Let us introduce the time operator: 

oo 

t= j t|t)(t|dt, (2.14) 

— oo 

having the following properties 

t|t) = t|t), (t|t = (t|t, t"|t) = t"|t). (2.15) 
In a way similar to the derivation of eq. (2. 13) we can get: 

(«l(-it)|u> = ^(«|u>- (2.16) 



From eq. (2.7) we obtain 



dt 

and more generally (rewriting eq. (2.13): 



i— (t|co) = co(t|co), (2.17) 



-i^(t|u) = (t|a)|u), (2.18) 

i.e. according to eqs. (2.17) and (2.18) one can consider the operator 

-i— = c5, (2.19) 
dt 

also as the angular frequency operator. The difference between the operators of eq. 
(2.11) and eq. (2.19) is that the operator co is an infinite dimensional matrix in the Dirac 
space, while the operator w acts in the space of functions. In the following we will deal 
mostly with operators in the Dirac space. 
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3. Continuous Spectral Analysis 

The spectral decomposition of a time dependent signal f(t) is given by its Fourier 
transform: 

+00 

FN^vfel e-'"f(t)dt, (3.1) 

—00 

while the inverse relation is: 

f(t)=^ j e'-F(co)dco. (3.2) 

—00 

Here we will utilize the Dirac bra-ket formalism [9], (again, in the rigged Hilbert space 
[18]), for spectral analysis. In practice the data are collected in finite samplings. In this 
section, for methodological reasons, we will confine ourselves to continuous signals. 

We will denote the signal function f(t) as a bra-ket (t|f). Thus the signal is 
considered to be an abstract ket-vector state |f) and can be represented either by (t|f) 
or by its Fourier transform (co| f ) : 

+00 

((o|f)=j ((o|t>(t|f>dt, (3.3) 



and the inverse transform is 



(t|f)= j (t|co)(co|f)dco, (3.4) 



where: 



F(co) = (co|f) , (3.5) 



and the normalization condition is: 

00 00 

(t|t')= j (t|co)(co|t')dco = ^ I e^^^^'^'ho)^ 8(t-t'), (3.7) 

— 00 — 00 

where is the delta function of Dirac. In the Dirac formalism the completeness relation 
is: 



i=j |t)(t|dt=j |co)(co|dco (3.8) 

— oo — oo 

where I is the identity operator. The relations (3.6) are solutions of the eigenvalue 
equation (2.17) 

3a. Windowing 

In practice the signal is detected during a finite time interval, say 2T. The spectral 
analysis in this time interval can be done in the following way: 

T 

F(co)T = (co|A|f)^ =^ j a(t)exp(-itco)f(t)dt (3.9) 

-T 

where a(t) is one of possible windowing functions [22-25], [29] and eq. (3.9) was 
obtained by using the filter operator 

T 

A = I a(t)|t)dt(t|. (3.9a) 

-T 

In this case the formalism can be replaced by decomposing the identity operator into: 

i = A + (i - A) = j a(t)| t)(t |dt + j [l - a(t)]| t)(t |dt , (3.8a) 

— oo — oo 

where a(t) is the windowing function having support only in the finite interval. For 
example for the rectangular window with a finite time interval 2T: 

a(t) = 1 , for -T < t < T ; (otherwise a(t) = 0.) (3. 10) 

One may employ also other windowing procedures of common use. Substituting eq. 
(3.8a) into eq(3.3) we obtain: 

+00 +00 

(co|f)=j a(t)(co|t)(t|f)dt+ j [l-a(t)]{co|t){t|f)dt 

— oo — oo 

+00 +00 +00 

= j a(t)(Q)|t)(t|f)dt+ j dco'j [l-a(t)](co|t)(t|co')(co'|f)dt 

— oo — oo — oo 

+00 

= (co|f)^+j K(co,co')(co'|f)dco', (3.11) 

— oo 

where 



+ oo 

K(co,co')= j [l-a(t)]{co|t)(t|co')dt (3.12) 

— oo 

is the kernel of the linear integral equation (3.11), and 

+CX5 

(co|f)^= j a(t)(co|t)(t|f)dt (3.13) 

— oo 

is the spectrum of the windowed signal (the measured quantity). Equation (3.11) can be 
recast in an operatorial form if we substitute: K(a),co') = (co|K|co') = (a)|l-A|a)'), and 
thus obtain from eq. (3.11): 

(co I f ) = (co I f )^ + I (co |k| co')(co' I f )dco' . (3.14) 

— oo 

Eq.(3.14) can be used for an approximate estimation of the error of the spectral 
analysis in the following way: 

(co|f)-(co|f)^- j (co|K|co')(co'|f)^dco', (3.15) 

and when possible (i.e. if eq. (3.14) has a solution), to get the exact spectrum form the 
windowed one. In the Appendix B we give an example of a continuous signal with 
windowing. 

3b. Filtering 

The filtering can be described with the help of the filter operator 

oo 

4)= j (p(co)|co)dco(co|, (3.16a) 

— oo 

which is an operator valued , continuous , linear functional. The filtered spectrum is: 

oo oo 

(co|4)|f) = (p(co)((o|F)- j (co|t)(t|4)|f)dt= j ((o|t)(t|f)^dt, (3.16b) 

— oo — oo 

where the function (p(co) modifies (filters) the spectrum and (t|f)jp = (t|4>|f) is the 

filtered signal. In a similar way to that derived in sec. 3a, one can derive the obvious 
integral equation which relates the signal to the filtered signal: 

oo 

t|f) = (t|<i)|f) + {t|(l-4))|f) = (t|<i)|f)+ j {t|(I-6)|t')(t'|f)dt'. (3.16c) 
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Let: 

oo 

(p(co)= j (co|t)(t|(p)dt, (3.17) 

— oo 

then, using eq. (3.17), the inverse transform of (3.16) is 

oo oo oo 

(t|f)^=j (t|(o)(p((o)(a)|f)d(0 = j j (t|(o)((o|t')(t'|(p)((o|f)doxit'. (3.18) 

—oo — oo — oo 

From eqs. (3.6) we have: 

(t|co}((o|t')=(t-t'|co), 
and eq. (3.18) becomes, after using the r.h.s of eq. (3.8): 

oo oo 

(t|f)^ = I (f |(p>(t-f|f>dt' = I (t-f|(p)(f|f)df (3.19) 

— oo — oo 

in which we recognize the convolution theorem. The inverse of eq. (3.17) is 

oo 

(t|(p)= j (p(co)(l/V27t)exp(itco)d(0. (3.20) 

— oo 

Let us consider for example: 

0, otherwise, 

than 



1 for coi <co<c09 and -a)9<(o<-0Di 
(p((0) = \ . ^ ^ ^ (3.21) 



— OO ^ 

If the frequencies of the signal are expected to be in the above mentioned region, i.e. 
bounded by eq. (3.21), then (t|f)^p ^ (t <i>|f) = (t|f) and eq. (3.22) becomes an identity 



[sin((t-t')a)2)-sin((t-t')a)i)](t-|j>|f) ^^,^ ^^^^^^ 



which is a particular case of eq. (A. 9) in Appendix A. Later we will relate this result to 
the sampling theorem. 



3c. Averaging 

One can further extend the Dirac formalism and obtain easily both well known, 
and as well, less known results. For example: 
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(f|f>=j (f|t)(t|f)dt= j (f|coXco|f>dco, 



(3.23) 



is obtained directly from eq. (3.8) and is the well known Parseval theorem. One can 
define average quantities in a manner similar to that done in quantum mechanics, e.g.: 



oo oo 

n \ \ I I / — oo — oo 



w = 



(f|f> (f|f> 

or work directly with the signal f(t) using eq. (2. 18): 



(f|f> 



(3.24a) 



j (f|t)(t|cb"|f)dt 1 (f|t)(-iA)"(t|f) 



(f|f> (f|f> 
In order to define an average signal we introduce the signal operator: 

oo 

f= j f(t)|t)(t|dt, f|t) = f(t)|t), 

— oo 

with the above definition the average will be 



(3.24b) 



(3.25) 



f = 



f|f|f 



An analog of mean kinetic energy will be 



K = 



dt 



t f 



> ^ 

/ -oo 


mm 


;t|f>dt 


) 


m 




oo 

J(fl'> 

— oo 


-m 
.lit . 


'(t|f)dt 


(flf> 



(3.26) 



(3.27) 



We can see that the introduction of the Dirac formalism in Fourier and signal analyses 
can lead to simplifications and to new definitions of averages. 

3d. The Autocorrelation and the Wigner Double Distribution Functions 

Let us consider the square of the absolute value of the signal S(t): 



t|S)|' =(S|t)(t|S)= j j (S|co)(co|t)(t|co')(co'|S)dcodco' 

— oo — oo 

= — f f (S|co)(co1S)exp(it(co'-co))dcodco' 



(3.31) 
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and substitute co = co' + x, co' = - x / 2 , we obtain: 



^ oo oo 

t|s)|'=— j j ^s|i^ + ^^(JiJ-^|s^exp(-itx))dMT, 



(3.32) 

eq. (3.32) introduces the point-splitting regularization which gives the Wigner function, 
which we derive next. Depending on the order of integration, we get two very important 
functions 



1 °° 

t\sf = — f R(x)exp(-itx))dx, 
271 •' 

— oo 

1 CO 

|(t|S)|'= — j W,(t,d)d^, 



(3.33) 
(3.34) 



where: 



S di^. 



— oo 

I oo 

(t, d) = — j (S I d + - I S) exp(- itx))dx , 



(3.35) 



(3.36) 



are the autocorrelation and the Wigner functions, respectively. Let us consider now the 
absolute value squared of the Fourier transform of the signal S(t): 



(3.37) 



|(co|S)| = (S|co)(co|S)= jj (S|t)(t|co)(co|t')(t'|S)dtdt' 

— oo — oo 

-t oo oo 

= — f f (S|t)(t1S)exp(ico(t-t'))dtdt', 

271 •' 

— oo— oo 

following the steps which led to eq. (3.34) we will get: 

|(co|S)|' = — [ W2(t,co)dt, 

271 

where: 



(3.38) 



oo 

(t,co) = — j (s|co + J^)(co-^ S^exp(-itx))dx. 



One can easily prove by using substitutions that: 

Wi(t,co) = W2(t,co), 



(3.39) 



(3.40) 
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therefore we have one double distribution function W(t,Q)), the Wigner distribution 
function, satisfying 

^ oo 

W(t, = ^ j (S I « + ^/2)(^ - 1 S) exp(- itcoO)dco' 



■t oo 

2n' \ ' 



(3.41) 



S)exp(- io>T))dx. 



and (from eqs. (3.34) (3.38)): 

.| oo .| oo 

|(co|S)|' =— f W(t,co)dt, |(t|sf = — [ W(t,co)dco. (3.42) 
111 111 

— oo — oo 

In ref (21) one can find a detailed exposition of the Wigner function and other double 
distribution functions. 

4. Non-orthogonal states in subspaces 



Let us note that by going to subspaces of time or frequency the orthonormality 
condition might not be satisfied any more. If we restrict ourselves to the windowed 
subspace, the rectangular window for example, we obtain: 



(a),T|a)',T)= j (a)|t>(t|a)')dt = | c'^^^-^^^dt= 



.1(0) -co)T_^-i((0--a))T _ sin(((o-coOT) 
27ri(a)-a)') 7t(a)-a)') 



(4.1) 



where 1®'^/ are the projected states into the above mentioned subspace, which we will 



denote by St. The orthogonality condition is satisfied only if 

co-co' = ±n7t/T, n=l,2,... 



(4.2) 



In order to write eq. (4.1) more explicitly one needs to introduce the projection operator 



1 |t>(t|dt. 



(4.3) 



which projects into the subspace St . The ket-vector states in this subspace now are: 

'|t), for|t|<T 
0, for Itl > T, 



t,T) = P^|t 



(4.4) 
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|co,T) = Pt|co)= j |t){t|co)dt= j |co')j {co'|t)(t|co)dtd(o' 



-T -~ -T 



17, , T -.tr. rv\ °? sinfTfco - CO')! 
— CO' e"("-"Wco'= dco'— ^ 

271 ' '-"^ TlfcO-CO') 

— oo —1 — oo \ / 



(4.5) 



CO 



Therefore eq. (4.1) should be written as follows: 



T 



co,T|co',T) = (co|Pt,Pt,|co') = 



> f (co|t)(t|a)-)dt= ^^"^^^-^^^^\ (4.6) 



The solutions of eq. (4.2), i.e. states which are orthogonal, are: 

«„,c =n7t/T + C, n = 0+l+2,..., (co„c,t|co^c'T) = 5„^T/7i, (4.7) 

where C is an overall constant and the normalization comes from eq. (4. 1). The states of 
eq. (4.7) form a complete orthogonal set, because we know that the functions 

{t|co„,o,T) = e'™"'°VV27r (4.8) 

form a complete orthogonal basis for Fourier series in the interval -T<t<T. Thus the 
ket- vector states |to,T) ^re an overcomplete basis in St . As we shall see later, the 
relation between the overcomplete basis |^'T) ^nd the complete basis of eq. (4.7) is the 

cornerstone of the sampling theorem. Let us first note that the projection operator is 

the identity operator in the subspace St : 

Pt|co,T) = PtPt|co) =.Pt|co) = |co,T) _ (4_9) 

Therefore we should expect the completeness relation, taking into account (4.7), to be: 

Pt=-L |«n,C'T)(co„^c'T|. (4.10) 

^ n=— 00 

Employing eqs. (4.6) and (4.7), the expansion of the overcomplete basis in terms of the 
complete one can be obtained in the following way: 

|co,T) = Pt|co,T) = ^ £ |co„,c,T)(co„,c,T|co,T) 



T„tl' ' 7i(co„c-co) ntl' ' T(co„c-co) 
from which the following sampling theorem, for a signal f, can be deduced: 
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- sin[T((0„p -co)l , 

this result was obtained under the assumption that the signal exists only for times t: 

(4.13) 



Before closing this section we will prove that the states \^^^/ form a basis which is 
overcomplete. We start with the completeness relation (3.8) and modify it according to: 

CO CO 

PT,=PT,iPT,= j PT,|co)(co|PT,dco= j |co,T)(co,T|d(0 (414) 

—00 —00 

Relation (4.14) is the proof that the states 1®'^) form a basis in the subspace St . 
Indeed any function with the support in St have a unique expansion: 

00 

u(t) = {t|u) = (t|PT|u)= j (t|co,T){co,T|u)dco_ (4,15) 

—00 

Next we multiply eq. (4.5) from both sides with and obtain: 

, °; sin[T(co - co')l , 
PJco,T)= co,T)= f dco'— L-^ ^ t^''T • (4 16) 

' ' ' _^ 7l(C0-C0 ) ' ' V • 7 

Eq. (4.16) shows that the states \^^^) are dependent between themselves, therefore the 

basis is overcomplete and according to ref (14) these states satisfy the conditions for 
being coherent states. 

5. Sampling theorem for band limited signals 

In section 4 the sampling theorem was derived for discrete frequencies for bounded 
time. In this section we derive the sampling theorem for time states in the case when the 
frequencies are restricted by the condition 

H^^B. (5.1) 

Condition (5.1) restricts the Dirac space to a subspace which we denote as Sb . This 
subspace is reached with the projection operator: 

COB 

Pb= 1 |(0>(co|d(0. (5.2) 

-COB 
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The ket- vector states in this subspace now are: 

. , , flco), for Icol < cOr 
co,B =P3CO = ' ' 

' ' ' ' 0, for |co| > cOg, 

B oo B 

|t,B) = PB|t)= j |co)(co|t)dco= j |t')j (t'|co)(co|t)dcodt' 



(5.3) 



-00 -B 



1 

1% 



V L ^t'-t) 



^. ^B(t--t)] ^^ 



(5.4) 



Let us check the scalar products 

COB 



(t,B|t',B)= j (t|a))(co|t')do) = 



-COB 



sin((t-t')C0g) 

7U(t-t') 



(5.5) 



the orthogonality condition is satisfied only if 

t-t' = ±n7t/C0B, n = l,2,... 

The solutions of eq. (5.6) are the times: 

t„D =n7t/C0B +D, ±n = 0,1,2,...^ 

where D is an overall constant. The corresponding states satisfy the orthogonality 
condition: 

and the identity operator for the complete set in Sb is the projection operator: 

PB=(^/«B)l|tn.D.B)(t„,I„B 

n=-oo 

The above results (eqs. 5.5-5.10) allow us to obtain 



(5.6) 



(5.8) 



(5.9) 



(5.10) 
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u(t) ^ (t|u) = (tlP^lu) = — £(t|t„^o,BXt„,D3 



CO 



u 



B n=-oo 



~ sin[(t-t„u)coB] 

= L ~uz, — \ — ^u(t„D), 

n=-~ ^^t tjj j3 jCOg 



(5.11) 



which for D=0 is the well known sampling theorem, and can be extended for D^tO. 
Eq. (5.11) can be recast into a more general form: 

TT ^ / M < ^ sin[(t-t„ n)©^] , 

(t,B| = (t|PB = ^ I(t3|t„,„B)(t„,o3| = I . T '\ t,,,B 



CO 



B n=-<x) 



(t-tn,DK 



(5.11a) 



or: 
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|t,B) = PB|t) = ;f i:|t,,,B)(t,o>B|t,B)= I °7 |tn,D3). (5.11b) 



®B n=-oo 



The inverse relation can be obtained in a similar way: 

" " sin[(t — t )co ] 

tn,D3) = PB|t„,D3)= j |t,B)(t,B|t„,D3)dt= j dt V_°' / |t3). (5.11C) 

Let us decompose the signal u(t) into its odd and even parts: 

u+(t) = i[u(t) + u(-t)], u_(t) = |[u(t)-u(-t)]^ (5.12) 

and let us denote 

tn=t„,o, (5.13) 

than, from eq. (5.11) the odd and even parts of the signal u(t) can be presented in the 
following way: 

^ sin[(t-tJ(OB] ^ sin[(t + t„)(OB] 

„=i (t-tJ(OB (t + tJCOB 

sin(tcOB) - sin[(t-tJcOB] ^ sin[(t + 1 JcOb] ^^'^"^^ 

u+ (t) = — (0) + ^u^ (K) + L V. — ^ 

In eq. (5.14) the signal is sampled only for non-negative times. This is an interesting 

possibility if the time parity is known and fixed. 

We can consider the sampling theorem from a different point of view if we 

substitute eqs. (5.8) and eq. (5.13) into eq. (5.11), we obtain: 

f (- 1)° sin(cDBt) 
u(t)= £ u(tj, (5.15) 



or: 



cOnt-nTi; 



u(,) ^ - J^„(,„). (5.16) 



sin(cOBt) n=— oo 

Thus eq. (5.16) is the pole expansion of its left hand side, and it is valid only if u(t) 
does not have singularities in t except at infinity. It has also to comply to the Mittag- 
Leffler 's theorem [26], 

We can still view the sampling theorem from other angle. In eqs. (5.6)-(5.10) we 
have different orthogonal bases characterized by the parameter D, let us call it B(D). 
Every point t can be reached by t=tn+D for some values of n and D. Therefore we can 
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present the sampling theorem from that point of view by expanding the base vectors of 
B(D) in terms of the base vectors of B(D=0). For this we need only to use eqs (5.11) . 
We obtain: 



tn,D 



\ Sin[(t„,D-tn,)«B] 



and the sampling theorem takes the form: 
t = t^ 



where now only discrete orthogonal sets are involved. Eq. (5.18) a consequence of 
isometric isomorphism between the bases labeled (k,D) and those labeled (k,0). 

It is interesting to note that all band limited functions, f (t) = (^IPbI f) , have to satisfy 
the identity given by eqs. (A. 10) and (A. 1 1) of Appendix A. Using eq. (5.5), we obtain 

CO 

(t|PB|f)= j (t|PB|t')(t'|P3|f)dt' 

7 sm(cOR(t-t )) , , , , 
i 7r(t-t') ^ I «l / 

The identity (5.19) was employed in refs. [27] and [28] in order to derive an expansion 
for band limited functions in a time window. 

6. Some consequences of the sampling theorem and some ambiguities 

First let us note, that in order to avoid some ambiguities, the condition (5.1) has to 
be changed to 

H<«B, (6.1) 
otherwise we can add to the eq. (5.11) any function which is zero at the interpolation 
points tn,D. Thus we can get an ambiguity of up to a function having zeros at t= tn,D. 
Such a function will have a periodicity of = ^'^b . The condition (6.1) will eliminate 
that possibility 
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A consequence of eq. (5.10) is the replacement of the Fourier transform eq. (3.3) by 
the sum 

(«>|f> = ^ t {<^K){K.\() = ^ t e-"-f(t„,o). (6.2) 
where we continue to use the normalization given by eq. (3.6) 

The sampling theorem is recovered from eq. (6.2) by the integration: 

COB 

(t|f)= j (t|co)(co|f)da)_ (6_4) 

-COB 

We shall derive, using eqs. (5.10) and (5.8), an important relation to be used later on 

OO -J CO 

5(a)-co') = (a)|a)') = ;f I {^K){K,oH = I. e'-("'-)^"B _ ^^ ^^ 

"Jb n=-~ ^"^B n=-~ 

The right hand side of eq. (6.5) is a periodic function of period , therefore we infer 
that it should repeat itself every period, thus we should have: 

-1 OO OO 

1 ^in^co/coB = £ 5(co - 2a)Bm) . (6.6) 

^^B n=-~ m=-~ 

The last relation is proved more rigorously in ref [3]. 

Caution has to be taken in using eq. (6.2) which has a periodicity of ^ = ^^b, 
in contradiction to eqs. (5.1) and (6.1). The restricted limits integration (6.4) bring us 
back to eqs. (5.1) and (6.1). Therefore eq. (6.2) can be used without ambiguities only 
under the condition (6.1). 

Let us consider an example: 

F(t) = (t|F) = e'"o\ a)o>a (real)^ (6.7) 

and let us check eq. (5.11) by first evaluating eq. (6.2) (with D=0), than eq. (6.4). 

/IT 00 00 

(co|F) = ^|^ e'("o-co)tn = 5(0)0 -a)-2m(0B), KI<«b. (6.8) 

where eq. (6.6) was used. Insertion of eq. (6.8) into eq. (6.4) recovers the correct result 
i.e. eq. (6.7). If instead of condition (6.1) the condition (5.1) will be used, than one is 
allowed to take = ^'^o and instead of eq. (6.8) we obtain: 
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(co|F> = ^^£ e^("o-co)t„^^2^ 5(0)3- CO -2mt03), 0)b=«o, (6.9) 
and after substituting in eq. (6.4) we obtain (the delta function is an even function): 

COB 

(t|F)= j e'™[5(o)-0)B) + 5(o) + 0)B)]dO)7^e^"o\ (6.10) 



-COB 

i.e. we do not recover the resuh (6.7). 

Now we shall derive another application of eq. (6.6). Substituting the Fourier 
transform (3.2) into eq. (6.2) we obtain: 

(6.11) 

= F(«-2cOBn), 

n=— 00 

which coincides with eq. (6.2) under the constriction (6.1). If no constraints are 
imposed on the frequencies, eq. (6. 1 1) and its complex conjugate are the manifestation 
of the aliasing phenomenon. 

The main ingredients leading to the sampling theorem are eqs. (5.5), (5.9) and (5.10). 
A consequence of eq. (5.9), using eq. (5.5), is 

00 

t„,D|tm,D)= j (t„,D|t)(t|t^,D)dt = 6^0)B/7i. (6.12) 



or 



sin[(t„D -1)0)3] sin[(t^D -t)o)B]^ 5 

dt = 



nm 



(6.13) 



L (tn,D-t)COB (t„,,D-t)tOB CObTI- 

7. Spectral analysis on the half-infinite time axis. 

In the sampling theorem and related equations of section 5 the time is running from 
minus to plus infinity. In practice we start to collect data at a definite time, which can 
be chosen to be t=0. Therefore half axis time analysis (t>0) is more related to 
experiment than whole time axis analysis. Moreover we deal with real signals. For real 
signals the Fourier transform eq. (3.1) will have the following property: 

F*(co)=F(-co). (7.1) 
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Let us define the even and odd parts of the signal F(t) of eq. (3.2) as: 

f,(t) = i[f(t) + f(-t)], f_(t) = l[f(t)-f(-t)], (7.2) 

than, using eqs. (3.1), (3.2), (7.1), (7.2) we obtain the sine and cosine Fourier 
transforms: 

f^(t) = -^j cos(cot)i[F(co) + F*(co)]dco= -j cos(Q)t) Re[F(co)]dco , (7.3) 
and in a similar way: 

[2"+'" 

f_(t) = - -j sin(cot)Im[F(co)]dco, (7.4) 







Re[F(co)] = J- f cos(cot)[f^(t)]dt, (7.5) 



+00 



lm[F((o)] = - -j sin(cot)[f_(t)]dt. (7.6) 
' ^ 

In terms of the Dirac formalism the sine and cosine transforms can be formulated as 
follows: the vector space with the basis 1 1) is decomposed into two orthogonal 

subspaces | t,P = l) (with projection P=l into the even functions oft) and | t,P = -l) 
(with projection P=-l into odd functions oft). Now, because of the symmetry of even 
and odd functions in t ^-t, all the information can be stored in t>0. Thus any function 
oft can be decomposed into its even and odd parts, as in eq. (7.2): 

f (t) ^ (t I f ) = [f, (t) + f_ (t)] ^ [(t, P = 1 1 f ) + (t, P = -1 1 f )] , (7.7) 

from which we deduce that: 

|t) = |t,P = l) + |t,P = -l). (7.8) 
Let us note from eqs. (3.6), (3.7), (7.8), and (7.9) that: 

(t, P = 1 1 CO, P = 1) = cos(a)t), (t,P = -l|a),P = -l) = sin(a)t), 
J 2% Jin 



(t,P = l|t',= P-l) = 0, (7.9) 
(t,P = l|t',P=l) = (t,P = -l|t',P = -l) = i5(t-t'), 

In a similar way we can define even and odd parity states of frequency. 
The identity operators will acquire the forms: 



(7.10) 
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oo oo 

i= j|t,P=l)(t,P = l|dt+ j|t,P = -l)(t,P = -l|dt 

— oo — oo 

oo oo 

i= j|co,P = l)(co,P = l|dco+ j|co,P = -l)(co,P = -l|dco 

— oo — oo 

For real signals, eqs. (7.3)-(7.6) will acquire the following form: 

oo 

(t,P = l|f) = 2j (t,P = l|co,P = l)(co,P = l|f)dco (7.3a) 



oo 

(t,P = -l|f) = 2j (t,P = -l|co,P = -l}((0,P = -l|f)dco, (7.4a) 



fco,P = l|f) = 2j (co,P = l|co,P = l)(t,P = l|f)dt, (7.5a) 



co,P = -l|f) = 2j (co,P = -l|co,P = -l)(t,P = -l|f)dt (7.6a) 



We see in eqs. (7.3)-(7.6), (7.3a)-(7.6a) that in the sine and cosine transforms the 
integration is for positive time, but the real signals must have a definite parity, i.e. they 
must be either even (P=l) or odd (P=-l) functions oft. Let us now consider the general 
case with no definite parity. 

As we have pointed out, in an experiment we measure the signal, starting to 
collect data at a definite time, which can be taken as t=0. Before that time our 
apparatus is idle. Let us consider for simplicity the two signals: 

t|Fi) = -^exp(ia)ot) = (t|a)o), 

0, fort<0, (7.11) 
exp(itOot), foryt>0. 



t F2 



, /27C 

and evaluate their Fourier transforms, from eq. (3.3): 

00 

co|Fi)= j ((jL)|t)(t|(jL)o)dt = 6(co-cOo), 



1 " 11 

(C0IF2) = — j exp[i(cOo -co)t]dt = + ^6(cOo - co), 

271 27i;i CO - cOq 



(7.12) 
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Thus in a positive time spectral analysis a single frequency signal can be represented 
by a pole in the frequency with a possibility of using rational approximations. Let us 
introduce the positive time projection operator 



t^ =1 |t)(t|dt, (7.13) 







and apply it on the kets and bras of our vector space. We will note that the frequency 
states will no longer be orthogonal: 

(co|t>') = j («KXt|»> = ^,^ + ^5(co-co'). (7.14) 

Let us denote 

OO .| oo 

A(co) = (co |t^ I F) = I (co I t)(t |t^ I F)dt = -j= \ exp(-icot)F(t)dt , (7. 15) 



V27r 



for the signal F(t). The inverse relation (for positive t) is: 

CO CXD 

F(t) = (t|t^|F)= j (t|t^|co)(co|t^|F)dco = ^= j exp(icot)A((o)dco . (7.16) 



From eq. (A. 10) of Appendix A and eq. (7. 12) we infer that 

A(co) = lf 4^dco', (7.17) 



711 CO -co 

— OO 

which is a Hilbert transform for the real and imaginary parts of A: 

^ 17 Im[A(co')] , , 

Re[A(a))] = — — !— ^ — ^dco , 

TZ CO'-CO 

— OO 

^ -17 Re[A(co')], , 

Im[A(o))] = — — — -do) . 

% ^ ca'-o) 



(7.18) 



In eqs. (7. 15-7. 18) fast enough decay as co^oo is required for the existence of the 
integrals. 

Let us consider the case of the frequency bounded by eq. (6.1) and consider the 
sampling theorem for positive time signals. As in the case of eq. (6.2) one can replace 
the transform (7. 15) by the sum: 
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A(co) = {co|t+|F) = — {^\K 



^n,D /Vn,D 



n-i-r)>n'- 



'B n=-~ 



CO t t 



t„,D>dt 



(7.19) 



n,D 



'B n+D>0 

which is not the same as the (similar to eq. (6.2)) expression: 
B(co) = — £ (co 

®B n+D>0 



F> = i^ E e-«-F(,,„„). (7.20) 

V^lOg n+D>0 



Let us consider as an example the second signal of eq. (7.11) for which eq. (7.20), with 
D=0, becomes: 



B(co) = 



1 



2co, 



,i(COo-tO)n7c/a)B _ 



1 



1 



(7.21) 



'B n=o 2c0b 1- exp(i(cOo - co)7i; / ©b) ' 

which differs from the second expression in eq. (7. 12), although both expressions have 
the same pole and residuum. We mention this because eq. (7.20) is the basis for one 
sided Z-transform [5-8]. Instead, expression (7.19) should be used. Expression (7.20) is 
a good approximation to (7.19) in the sense that it has the same poles and residua. In 
eq. (7.19) there appear new coefficients: 



rcc , , w \ 1 pec sin((t-t„D)(OB) 
C(»,t.,„) = £ (0)|t)(t|t„,>dt = -^|„ e-^^^-^dt, 

this expression can be brought, with the help of eqs. (4.6) and (7. 14) to the form: 



(7.22) 



, X 1 "f exp(-ico't„n) 



2%_i^ 27ti(co'-co) 



(7.23) 



8. Summary and conclusions 



We have demonstrated in this paper that the Dirac representation theory can be 
effectively adjusted and applied to signal theory. The advantages of the Dirac 
representation theory is its clarity, transparency and generality. It is suitable to be used 
with generalized functions and general transforms. The use of the identity operator and 
projection operators simplify greatly the calculations and derivations. Using these 
operators we introduced in Appendix A, a large class of functions, which were called 
the incomplete delta functions, as they have similar properties to the Dirac delta 
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function and converge to it in a limiting process. They have been used throughout the 
paper in order to simpHty derivations. The case of wavelets deviates somewhat from 
the standard procedure, and is collected in Appendix C. 

The scalar products (brackets) of Dirac are written in a way which does not depend 
on the choice of the coordinates. In the process of representing the signal theory with 
the brackets we found interesting features of scalar products which were related to the 
change of validity domains of the conjugate coordinates (in signal theory they are the 
time and frequency, in quantum mechanics the coordinates and momenta). An 
orthogonal basis which spanned the time space ceased to be orthogonal after 
constraining the domain of frequencies (bandwidth) and becomes an overcomplete. But 
in this case there exist a sub-bases of discrete times which are orthogonal and 
complete. The completeness of the two bases, one for all frequencies (the overcomplete 
non-orthogonal basis), the other in the bandwidth, eq. (4.10), allows the derivation of 
the sampling theorem relating signals in the discrete time (complete) basis to signals 
for all times (with the overcomplete basis - eq. (4. 1 1)). Thus the relation between the 
overcomplete bases and a complete one is the essence of the sampling theorem. 
Moreover eq. (4.10) allows us to replace the continuous time integration by a discrete 
summation in all time integrations, this we could do in the Fourier transform, eq. (6.2), 
in the wavelets transform (Appendix C) and in the case of signals existing for positive 
time - eq. (7.19) . 

The idea of viewing the sampling theorem as the relation between an overcomplete 
basis and a complete one can be extended to quantum mechanics. We are in the process 
of completing a work in this direction. 

Appendix A 

The incomplete delta functions and the delta function of Dirac 

Let us consider two complete continuous bases: 

(x| x') = 5(x — x'), (y I y') = S(y - y') , — oo<y<oo; — oo < x < oo , (A. 1) 

linked together by the Fourier transform (eq. (3.6)) : 
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I y) = exp(ixy) . (A.2) 
V27r 



where 5(x)is the deha function of Dirac. The completeness relation is given by 

oo oo 

1=1 |x>(x|dx= j |y>(y|dy, (A.3) 

— oo — oo 

where I is the identity operator. Let us consider the following set of functions: 

5(x-x',a)= 1 (x|y)(y|x')dy 

—a 

= ±f e">dy =i!i^^fe^. a>0 (A.4) 
271 •' Tt ax 

— 3. 

from which the Dirac delta function 6(x) can be obtained in the limiting process: 

6(x)= lim6(x,a). (A.5) 

a— »cio 

The Dirac delta function (distribution, or generalized function) has the following 
properties: 

X 1- s / ^ I^' for x^O, 
o(x) = limo3(x) = < 

[oo, for X = 0. 

oo oo 

f 6(x)dx= lim f 5^(x)dx= 1, (A.6) 

— oo — oo 

oo oo 

f f(x)5(x-Xo)dx= lim f f (x)5^(x- Xo)dx = f (Xq). 

— oo — oo 

One can generalize the above procedure to a larger class of functions which we will 
call "incomplete delta functions". Let P be a projection operator, such that by changing 
its parameters, in the limiting process, it converges to the identity operator. An example 
is eq. (A.4) , in which the following projection operator was used: 

a 

P(a) = j |y)(y|dy; lim,^„P(a) = 1. (A.7) 

—a 

Let us define the incomplete delta function as 

5(x-x',P) = ((x|P)(P|x')), such that 5(x- x',P ^ I) = 5(x- x'), (A. 8) 



27 

i.e. the states denoted by x are projected to a subspace with the projection operator P. 
The projection operator satisfies (see eqs. (2. 1 la) and (2.23)): P=P and P=P', which 
can be employed in order to derive general properties of the incomplete delta functions. 
We can prove that similar equations to eq. (A. 6) are satisfied for the incomplete delta 
functions. Inserting the identity operator of eq. (A.3) into eq. (A. 8) we obtain: 



(A.9) 



5(x-x',P) = ((x|P) 1 (P|x'))= j (x|PP|x")(x"|PP|x")dx" 

— oo 

CO 

= j 5(x - x", P)5(x" - x', P)dx". 

— oo 

Let f(x) be a function with a support in the projected subspace, than: 

CO 

f(x) = (x|f) = (x|P|f) = (x|PlP|f)= j (x|PP|x')(x'|P|f)dx' 

(A. 10) 

oo 

= I 5(x-x',P)f(x')dx'. 

— oo 

For non restricted functions still the following identity holds: 

oo 

(x|P|f)= j 5(x-x',P)(x'|P|f)dx'. (A. 11) 

— oo 

The identities (A. 10) may have applications in testing hypotheses that a given function 
has a support in a restricted subspace. 
One can note from eq. (A.4) that: 

5(0,a) = a/7i. (A. 12) 

We shall prove the important property of the incomplete delta function (A.4): 

5(bx,a) = ^5(x,a|b|), (A. 13) 

|b| 

where b is a real number. For b>0 we have: 

, a , ab , 

5(bx, a) = — f e"'^dy = — f e"Mz = - 5(x, ab) . (A. 14) 

2%!^ 2%h_{^ b 

For b<0 we obtain: 

1 a 1 ab 1 1 

5(bx,a) = — f e"'^dy = — f e"Mz = f e"Mz = ^6(x,a|b|), 

which, together with eq. (A. 14), proves eq. (A. 13). From eq. (A. 13) we can get: 
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5(bx) = ^6(x). (A. 15) 

|b| 

Eq. (A. 15) can be generalized to include functions f(x) having zeros at x=Xn: 

n |f(xj| 

where near the zeros: 

f(x)»r(x„)(x-x„). (A.17) 
Appendix B. 

Analysis of respirator's cycle. 

The working cycle of period P of the respirator in the operation room of the 
Gynecology Department, Soroka Medical Center, Beer-Sheva, Israel, was 
approximately as follows: 

\/r t^ - J A sin(27i(p ^t) for < t < t ^ (pumping), . . 

^^^^"\0, for ti<t<P (idle), ^^'^^ 

where V was the volume of the pumped air, A is a constant volume, ty= 3 sec, 
(Pj=(l/6)Hz and the period of P was changing with the anesthetic procedure. Most of 
the time it was P=5sec, within the range 3.3sec<P<16sec. 

For practical calculation purposes it will be useful to translate t^t + P/ 2by half a 
period to obtain V(t) in a symmetric form for -P/2 < t < P/2 : 

vm- J^^^^^^^^il^~^2|) for t2<|t| (pumping), 
^'^^^"[O, for -t2<t<t2 (idle), ^ 

where 

t2=(P-t,)/2 . 

The above function is periodic, with a period P, and is an even function oft. 
Therefore a more restricted basis than (3.6) may be more suitable, namely the basis 
of even functions with period P: 

(t|n,P) = cos(2n7it/P), n=0,l,2,... , (B.3) 
and V(t) = (t I V) can be expanded in a following Fourier expansion: 
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V(t) = (t|V) = £(t|n,P)(n,P|v; 



n=0 



= V+£^nCOs(2n7tt/P), 



(B.4) 



n=l 



where (t|n,P) are periodic symmetric functions oft, orthonormal in the interval 
[-P/2, P/2] and 



^ p/2 

a„ = - j V(t) cos(2n7tt / P)dt 



t2 

2A 

TIP 



cos27r[(f„ -(pi)t + (pit2] cos27i[(f„ +cpi)t-q>it2] 



fn-9l 



f„+9i 



t=P/2 



t=t2 



(B.5) 



and 



(B.6) 



f„=n/P. 

Here fj=l/P is the basic frequency and f^ (n>l) are the higher harmonics. 
In practice measurements are taken for finite intervals. In this case the spectral 
analysis should be modified according to the procedure outlined in the main text 
eqs. (3.9) and (3.10). With the finite interval of duration 2T the spectral analysis of 
eq. (B.3) (with F(t)=V(t) in eq. (3.9) and remembering that CO = 27rf ) becomes: 



_ ao sin(coT) - J sin[(co - co JT] sin[(co + co JT] 
^v"J;t ~ 2j 1 ' 

CO n^i [ CO-COjj CO + COjj 

where CO = 27rf = 27rn/P. 

n n 

Appendix C 
Dirac Formalism as a framework for wavelets 



(B.7) 



We can show that the Dirac formalism can be applied to the wavelets transforms 
[29] of the signal f(t), which are defined in the following way: 

^t-b^ 



(w^f) (b,a) =|ar^ I mr 



dt 



V a 



(C.l) 
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Above in eq. (C. 1) \|/ is a square integrable function, the "basic wavelet". In order to 
apply the Dirac formalism we have to assume that the function \|/ is generated by an 
operator W which is represented in eq. (C. 1) by its matrix elements: 

(t|W|b,a) = |ar2 \|/f— Y (C.2) 
y a J 

Than eq. (C. 1), the wavelet transform of f(t), can be represented as: 

CO 

(b,a|W|f)= J (b,a|W|t)(t|f)dt. (C.3) 

— oo 

The discrete wavelet transform of Daubechies [30] can be treated in a similar way. Let 
us denote by Wd the operator generating the wavelets, with matrix elements : 

(t|Wo| j,k) ^ i|/j,,(t) = 2j^^i|/(2Jt-k) . (C.4) 
They form an orthonormal basis: 

J (l,m|Wi,|t)(t|Wo|j,k)dt =5j,-5,,^ , (C.5) 

— oo 

with the identity operator: 



1= L W^|j,kXj,k|Wi, . (C.6) 

j,k=- ~ 

The signal f(t) can be expanded using eq. (C.6) 

f(t)^(t|f>= L (t|Wi,|j,k)(j,k|Wi,|f) , (C.7) 

j,k=- °° 

where : 

(j.k|Wj,|f)= j (j,k|Wj,|t)(t|f)dt . (C.8) 

— oo 

Other discrete wavelet transforms [31] can go along an identical path. We shall now 
look for operators which can generate the wavelet transformation of equation (C.2). Let 
us assume that the operator V(a,b) has the following property: 

V(a,b)f(t) = -)=ff— 1 (C.9) 
Va V a ; 

Such an operator can be formed from the following operators [21]: 

V(a,b)=e-'('"''^^e-**, (C.IO) 
where C is the scale operator and C> the frequency operator defined by: 
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C = 



dt dt 



i dt 



The result (C.9) is obtained by translation: 

e-''''*f(t) = f(t-b), 

and by compression [21]: 



1 JO 



■JsL 



(C.ll) 



(C.12) 



(C.13) 



The transformation V of eq. (C.IO) is unitary and thus the normalization is preserved 
under its action. Using results of ref (10) we can represent the completeness relation 
and the inverse transform in the following way. The completeness relation can be 
represented by the identity operator: 

l = ^jj^|b,a)(b,a|^^, (C.14) 



where the normalization constant is defined by: 



|¥(co)|' 
Icol 



-dco <oo, 



(C.15) 



and vi'(co) is the Fourier transform of \|/(t) of eq. (C.2). Using eq. (C. 14) we can get the 
inverse transform of eq. (C.3): 

db da 



f(t) = (t|f) = -^jj (t|^|b,a)(b,a|^|f) 



In a similar way the wavelet transform (2.19) can be modified to: 

oo 

(b,a|^|f) = (71/0)3) Li^'^l"^ 

n=- 00 

and the wavelet transform (2.24) to: 

00 

( j, k |Da| f ) = (71 / CO B ) 2^ ( j, k |Da t „ D )(t n,D 



(C.16) 



tn,D/\tn,D 



(5.3a) 



(5.3b) 



32 



References: 

1. A. Gersten, Dirac's Representation Theory as a Framework for Signal Theory. L 

Discrete Finite Signals, Ben-Gurion University preprint, denoted as paper [I]. 

2. C.E. Shannon, Bell Syst. Tech. J. 27, 379-423, ibid 623-656, 1948 

3. J. Von Neumann, "Mathematical Foundations of Quantum Mechanics", 

Princeton University Press, Princeton NJ 1956 (translation from the German 1932 
edition) 

4. G.G. Emch, "Algebraic Methods in Statistical Mechanics and Quantum Field 

Theory". John Wiley, New York 1972 

5. S.L. Braunstein, CM. Caves and G.J. Milbum,y4««. Phys. (N.Y.) 247, 

135-173, 1996 

6. L. Schwartz, "Theorie des Distributions, vols. 1-2, Hermann, Paris 1950-1951 
7. 1.M. Gel'fand and G.E. Shilov, "Generalized functions", vols. 1-2, Academic 

Press, New York 1964, 1968 (Translation from Russian). 

8. M.J. Lighthill, "Introduction to Fourier Analysis and Generalised Functions", 
University Press, Cambridge 1958 

9. P. A.M. Dirac, "The Principles of Quantum Mechanics", 4th ed. Clarendon Press, 
Oxford 1958, (the first edition appeared in 1930). 

10. 1.M. Gel'fand and G.E. Shilov, "Generalized functions", vol. 3, Academic 

Press, New York 1967 
11. 1.M. Gel'fand andN.Ya. Vilenkin, "Generalized Functions", vol. 4, 

Academic Press, New York, 1964 
12. 1.M. Gel'fand, M.I. Graev andN.Ya. Vilenkin, "Generalized Functions", 

vol.5. Academic Press, New York, 1966 

13. K. Maurin, "General Eigenfunction Expansions and Unitary Representations 

of Topological Groups", Polish Scientific Publishers, Warsaw 1968 

14. N.N. Bogolubov, A.A. Logunov and I.T. Todorov, "Introduction to 
Axiomatic Quantum Field Theory", Benjamin, London 1975. 

15. J.E. Roberts, J. Math. Phys. 7, 1097-1104, 1966 

16. J.P. Antoine J. Math. Phys. 10, 53-69, ibid. 2276-2290, 1969 



33 

17. A. Bohm and M. Gadela, "Dirac Kets, Gamov Vectors and Gel'fand Triplets. " 
Springer, New York 1989 

18. S.J.L. van Eijndhoven and J. de Graaf, "A Mathematical Introduction to Dirac's 
Formalism", North Holland, Amsterdam 1986 

19. O. Pelc and L.P. Horwitz, J. Math. Phys. 38, 115-138, 1997. 

20. A. Friedman, Int. J. Theor. Phys. 34, 307-338, 1994. 

21. L. Cohen, "Time-Frequency Analysis", Prentice Hall, Englewood Cliffs, 
New Jersey 1995 

22. L.E. Franks, "Signal Theory", Prentice Hall, Englewood Cliffs, 
New Jersey 1969 

23. J.G. Proakis and D.G. Manolakis, "Digital signal processing: principles, 
algorithms, and applications", 2"'^ ed., Macmillan, New York 1992 

24. A.V. Oppenheim and R.W. Schafer, "Discrete-Time Signal Processing", 
Prentice Hall, Englewood Cliffs, New Jersey 1989 

25. M.B. Priestly, "Spectral Analysis and Time Series", Academic Press, 
New York 1981 

26. K. Knopp, "Theory of Functions", Dover, New York 1947 

27. D. Slepian and H.O. Pollak, Bell Syst. Tech J. 40, 43-64, 1961 

28. D. Slepian, Proc. IEEE 64, 292-300, 1976 

29. C.K. Chui, "An Introduction to Wavelets", Academic Press, New York 1992 
30. 1. Daubechies, "Ten Lectures on Wavelets", SIAM, Philadelphia PA 1992 
31. 1. Meyer, "Wavelets: Algorithms and AppHcations, SIAM, Philadelphia 1993 



